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ABSTRACT 


It has been observed that flares with greater than 10 Mev gamma- r/ay 
emission are concentrated around the solar limb with a dispersion of 10 to 
20 degrees. It is shown that the bremsstrahlung by relativistic electrons 
is responsible for such gamma-rays and that the expected relativistic 
beaming cannot explain this dispersion. It is argued that thi3 dispersion 
is predominately a reflection of the pitch angle distribution of the 
electrons. Then it is shown that this requires a small variation of the 
magnetic field from the point where the electrons are injected to the 
photosphere and a nearly isotropic (in the downward direction) pitch angle 
distribution at the injection. The influence of other effects on the 
observed distribution is also briefly discussed. 


I. INTRODUCTION 


The gamma-ray detectors (Chupp ot al 1901) on board the Solar Maximum 
Mission (SMM) have observed gamma-raya of up to 30 Mev from a few 3olar 
flares. The bulk of the gamma-raya (from about one to seven Mev) are 
produced by energetic protons (Ramaty et al 1983). The residual continuum 
at these and higher (10 to 30 Mev) energies must be generated by some other 
processes. A recent analysis of all flares with detectable signals at 
photon energies greater than 10 Mev by Rieger et al (1984) shows a strong 
concentration of bursts toward the solar limb. The histogram in figure (1) 
summarizes this result, where we have included four additional events from 
a more recent analysis (Forrest et al 1984). 

This distribution is a reflection of not only the directivity of the 

emission but also the intrinsic distribution of the flux (or the luminosity 

function) at energies greater than 10 Mev. Furthermore, because the sample 

is chosen from a larger sample that is limited to flares with fluxes 

greater than certain threshold value at energies of about 300 keV, the 

directivity of emission and the luminosity function at this lower energy is 

also a contributing factor. We shall first ignore the effects of the 

directivity at 300 keV (which is expected and has been observed to be 

smaller than that at 10 Mev, Vestrand et al 1984) and the luminosity 

functions. We shall return to these effects at the end of the paper. Then 

v 

according to the above histogram, the primary factor, namely the emission 
process, must be highly directional. The best candidate for this process 
is Coulomb bremsstrahlung by non- thermal electrons (contribution from decay 
of pion apparently is negligible, Ramaty et al 1983), which are known to be 


responsible for production of hard x-ray3, Then, the observed photons aro 
producod by electrons of energy greater than 10 Mev Lorentz factor 
y > 20) so that the relativistic beaming i3 very strong. 

The directivity of brerasstrahlung radiation in connection with solar 
flares have been discussed by various authors (Elwert and Haug 1971, Brown 
1972, Petrosian 1973). These studies were limited to photon energies below 
one Mev and ignored the curvature and convergence of the magnetic field 
lines, and the dispersion in pitch angle of the electron beams was treated 
in an approximate manner. It was clear, however, that the brerasstrahlung 
directivity i3 already significant at energies greater than a few tens of 
keV and increases rapidly toward higher energies. 

, More recent work (Leach and Petrosian 1982 and 1983, hereafter 
referred to as LPI and LPI1), where the effects of the curvature and 
convergence of field lines and the diffusion in pitch angle of electrons by 
the elastic Coulomb collision are included, shows that the low energy 
x-rays will not be highly directional (cf, also Leach 198-4 ) - This is 
because for the non-relatlvistlc electrons, which are responsible for hard 
x-rays, the energy loss rate and the diffusion rate in pitch angle are 
comparable so that even electrons with initial pitch angle equal to zero 
acquire a large mean pitch angle before they have lost a substantial 
fraction of their energy. In addition, the x-rays emitted by the lower 
energy electrons with their shorter mean free paths may come from the\ top 
as well as the lower portions of curved flare loops. 

Such isotropization of electron distribution and the consequent loss 
of directionality of the emitted photons becomes less and less pronounced 
at higher and higher energies because i) the diffusion rate in pitch angle 
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of electrons becomes smaller than the energy loss rate at higher energies, 
and 11) higher energy eleatrons, because of their longer mean free paths, 
penetrate deeper into chromospheric (or photospheric) layers where the 
magnetic field is presumably vertical and constant over the short radiating 
path. We will, therefore, first consider the expected bremsstrahlung 
radiation from relativistic electrons injected in a uniform magnetic field 
region ( § II). Then, in § III, we discuss the implication of the 

observation for the distribution of the accelerated electrons and the 

parameters of the flare plasma. The effects of other characteristics on 
the distribution are discussed in I IV and a summary is presented in § V. 

II. BREMSSTRAHLUNG RADIATION FROM RELATIVISTIC PARTICLES 

As Is well known, the bremsstrahlung radiation from relativistic 
particles is highly directional and is limited to within an angle 

0 1/Y th e direction of the emitting electrons. The cross section 

is a complicated function of photon energy k (which we will express in 
units of electron rest mass energy), electron energy , and the angle 0 
(see e.g. Koch and Motz 1959), but at small angles the angular dependence 
has the simple form (1 + 0 2 y 2 )“ 2 . At larger angles (0 order of unity) the 

cross section flattens out but at a level of y -4 below that of small 

angles. For the expected steep electron energy spectra, gamma-rays of 
energy k are emitted mainly by electrons of energy y- 1 slightly larger 
than k. Consequently, the photon distribution is also expected to be of 
the form (1 + k 2 0 2 )" 2 . As shown in figure 1, this form is much 

narrower than the observed distribution of the gamma-ray flares, indicating 
that the observed distribution is a reflection of the pitch angle 

distribution of the electrons (and perhaps other factors), and that to 
first order we can approximate the cross section by a delta function. For 


Page 4 


relativistic electrons of pitch angle a (or y = cos a) with respect to 
the magnetic field in a hydrogen plasma, this approximation gives (see eq. 
3BS(6) of Koch and Mots 1959) 

d 2 ff(k,0) D (32TTr 2 a /3)g(k,Y)S(cos0-y)dcos0dk/k , 
o t 

(1) 

g(k,y) - [1 + k/ Y f (3/4)(k/Y) 2 ]ln[2 Y (Y- k)/ke*] , 
where 0 now is the angle between photon momentum and magnetic field, a p 

— 1 3 

is the fine structure constant and r Q = 2.8 x 10 cm is the classical 

radius of the electron. Note that this equation is valid over a limited 
range of y/k ^ 1 . At large values of y/k one must use the screened form 
of the cross section, equation 3BS(6) of Koch and Motz. This however, 
alters only the logarithmic term and will not be important in our analysis. 


Now if F(Y,P,s) is the distribution of the electrons so that Fdydyds 
is the flux (electrons/sec, i. e. flux integrated over the cross section of 
the magnetic tube) in the energy range dy , pitch angle cosine range dy 
and the length interval ds along the field lines, then the flux of photons 
J(k,0)dkdcos 0 (photon/sec) is 


J(k, 0) = (327rr Q 2 a F /3k) 


00 , oo 

g(Y,k)dy F(Y,cos0,s)n(s)ds , 
1+lc 0 


( 2 ) 


where n(s) is the density of the ambient hydrogen atoms or protons. 
Following LPI, we define a dimensionless column depth x so that this can 
be written as 


J(k, 0) = (8a p /3klnA) 


03 

g(k,Y)G(Y,cos0)dY , 

. 

1+k 


G(Y,M) = 


F(Y>U,T)dx, dx .= 4irr lnAn(s) 


( 3 ) 


o 
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which shews that the angular distribution of the gamma-rays is determined 
by the integrated (over the whole emitting region) pitch angle distribution 
of the electrons. 

For a uniform and vertical field directed toward the center of the 
sun, the angle 0 varies from tt to tt/2 for flares from the center to the 
limb of the sun. This means that the relevant range of the pitch angles is 
n/2 < a < tt (or 0 < y < -1). F< * a more realistic case of electrons 
injected in a closed magnetic loop, the angle 0 in the right hand side of 
equation (3) must be replaced by a function £(0,*.), which depends on the 
shape and orientation of the loop as well as its location on the sun (of. 
e.g. Petrosian 1982), For example, for a semi-circular loop of radius R 
located in the plane of the solar equator £ = 0 + s/R - tt/2 where s is the 
distance measured from the top of the loop. However, as we mentioned above 
(see also below), it is unlikely that any significant number of gamma-ray 
photons can be produced from the upper portions of the loop where such 
complications must be taken into aacount. It should be emphasized, 
however, that the distribution of gamma-ray flares aoross the solar desk 
tells us only about the electron distribution with pitoh angles a > tt/2. 


III. TRANSPORT OF RELATIVISTIC ELECTRONS 

The kinetic equations for a steady state transport of electrons were 
developed in LPI and were applied to the x-ray emission from 
non- relativistic electrons in LPII (cf. also Leach 1984). We now extend 
the application of these equations into the relativistic regime. Let us 
first consider all of the relevant interactions of relativistic electrons 
Injected in a cold plasma of density n and of magnetic field of strength B. 
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These interactions are Coulomb collisions, brems3trahlung, synchrotron and 
inverse Compton. Table 1 shows the rate of energy loss and the 
characteristic photon energy for each pro;ess. Nearly, for production of 
>10 Mev photons, bremsstrahlung is the only relevant emission process 
because the others will require electrons of unreasonably high energies. 
However, the Coulomb collisions (at high n and low B) and the synohrotron 
radiation (at low n and high B) are the dominant energy 1038 mechanisms for 
relativistic electrons and must be taken into account. The value of the 
ratio 

R = VY C = 1 - 3(B/10 ° gauss) 2 (Y/20) 2 (20/lnA)(10 10 cm" 3 /n)(l-ii 2 ) (4) 

determines which of these processes will be the dominant loss process. For 

“10 Mev electrons at typical coronal values of B and n, H % 1, but at 

cor 

photoapheric values (n£ 10 16 em -3 , B = 1000 gauss) R phot £ 10“ 4 . 

The inverse Compton process can be ignored even as a loss process 
because, for normal solar soft photon density £. = Lq/HttRq c = 2 erg/cnf 3 

it is negligible as compared to synchrotron losses when the field strength 
is high, (Y g /Y IC ) = (B/7gauss) 2 , and it clearly is negligible as compared 
to Coulomb losses when the plasma density n is high. The inverse Compton 
process could become important only if the soft photon energy density £y is 
enhanced by three or four orders of magnitudes above the normal density. 

We therefore consider the effects of synchrotron and collisional 
losses on the distribution of relativistic electrons accelerated in a 
plasma of density n and field strength B. For this and especially for 
discussion of angular distribution of any emission (bremsstrahlung here) 
from such electrons, we need the changes these processes inflict on the 
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pitch angle of the electrons. These changes are derived in the appendix 
and the result presented as a rate of change of pitch angle cosine y /y in 
Table 1 . As evident, for relativistic particles, these rates are smaller 
than the corresponding energy loss rates. This implies that in a uniform 
field such particles spiral along ftSie field lines at a constant pitch angle 
till they become non-relativistic. The distance along the field traveled 
then is 


s = 
0 


ycgdy/ (Y s +Y c ) S PX 10 14 cm(10 10 cm" 3 /n) (Y/lnA) tan" 1 (it*) /R* 5 , 


(5) 


J 


Y 


which is much longer that the length of the coronal portion of the loop 
(with typical size of 10 10 cm, B s 100 gauss and n = 1O i0 cm -3 ) except for 
electrons with nearly 90 degree pitch angles (y < 10 -t| ). This means that 
relativistic electrons will penetrate through a column depth of 
N = jnds s 10 24 cq- 3 (Y/lnA) and lose most of their energy through Coulomb 
collisions below the chromosphere at densities n > 10 16 cm 3 . The yield of 
synchrotron and bremsstrahlung radiations will then be 


Y s - ft/3 - 3 x 10~ 5 (B/10 3 gauss) 2 (10 16 cm“ 3 /n)(y/20) 2 (20/lnA)(l -y 2 ) 
Y b e 5 x 10" 4 (Y/lnA) (lny/3) 


( 6 ) 


which are larger than the hard x-ray yield of about 10~ s (of. e. g. 
Petrosian 1973) of non-relativistic electrons. The collisional loss time 
is of order of a millisecond, which is an order of magnitude shorter than 
the time for traversing a loop of length 10 I0 cm and is much shorter than 
the duration of a flare ("10 sec). This means that the acceleration of 
relativistic electrons (like that of the x-ray producing non-relativistic 
ones) must continue throughout the duration of the gamma-ray bursts. 
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If the accelerated electrons are injected in the corona at the top of 

a loop, they will encounter converging field lines and their pitch angles 

will change according to the adiabatic invariance relation 

(1 - y z )/B s constant. Then electrons with pitch angle in the range 

7T/2 > a > sin (B„/B will be trapped above the transition region, 

u cr • 

where B q and B fcr are the strengths of the field at the top of the loop and 
at the transition region, respectively. 

Electrons with initial pitch angle smaller than sin" 1 (B q /B )^ will 
penetrate below the chromosphere, where the density soale height is much 
smaller than the magnetio field variation scale and radiate as in the 
uniform field case described above. The trapped electrons, on the other 
hand, will bounce back and forth till they have lost most of their energy 
(and have acquired non-relativistic energies) via synchrotron radiation (if 
R > 1) or Coulomb collisions (if R < 1). This process will last for a 
period of 

t Q - S q /c “ 3 x 10^s(10^cm ^/n) (y/lnA) tan ^(R^)/R^ (7) 


with synchrotron and bremsstrahlung yields of the order of 
, , t ( R - 3R 2 /5 + .. . R « 1, 

Y -l-tan-V)^ 

S 11- TT/2R'“ + . . . R » 1 , 


Y b a 5 x 10“ 4 (y/lnA)(lny/3) [ln(l + R)/R] 


( 8 ) 


This is different than the situation for non-relativistic (x-ray producing) 
electrons whose pitch angle changes on a time scale comparable to their 
energy loss time scale. This allows them to penetrate below the transition 
region and lose energy more quickly and increase their bremsstrahlung yield 
relative to their synchrotron yield. 
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The long life time and the low bremsstrahlung yield of the 

relativistic electrons dictates against a trap model unless the density in 
the coronal part of che loop is >10 I2 cm or R > 100 (e.g. n = 10 lo cm , B 

> 1000 gauss). The first condition seems unlikely and the seoond one 
implies even a lower bremsstrahlung yield (Y " 1 , Y b < 10“ s ). This would 

mean a total gamma-ray flux that is 10 s times smaller than the total 
microwave and submillimeter fluxes which is contrary to observations. We, 
therefore, conclude that the fraction of relativistic electrons with 
initial pitch angles a > sin -1 (B /B^)”* mu3t be small. This is 

satisfied if either the field is nearly uniform C B q B tr ) and/or the 
accelerated electrons (injeoted at the top of the loop) have small pitch 
angles. 

For the purpose of illustration, let us assume that the injected pitoh 

angle distribution is gaussian F^CYtjj) « exp{-(1 -y 2 )/a 2 }. Then the 

above discussion implies that the cases with a 2 » B /B„ can be ruled out 

o o cr 

because they will give too many trapped electrons. 

If a 2 « which means a highly beamed distribution of Injected 

electrons, then very few electrons are trapped. However, in this case, 
even below the transition region the electrons will be predominantly beamed 
downward and most of their bremsstrahlung radiation will be directed into 
the photosphere and absorbed there (as the albedo of the photosphere for 
>10 Mev photons is negligible). 

The analytic expression derived in LPI are applicable in thi3 case. 
As there is negligible Interaction above the transition region, the pitch 
angle distribution at the transition region will be a broader gaussian 
F t r(y»y) * exp{-(1 - y 2 )B Q /a 2 B tr ) than the injected one (of. LPI eq. 
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7). Below the transition region the fiold is uniform, the density is high 
(so that R « 1 and synchrotron losses are negligible), and the small pitoh 
angle solution in LPI can be easily extended to give (of. equation A, 9 in 
the Appendix) 

F(y,u,t) “ 2F Q (Y+T)exp {-a 2 /(a 2 +£)}/(a 2 +5), C - y~ 1 t/(y + t) . (9) 

This shows that throughout the whole loop the dispersion in the pitoh angle 

is less than a + 2/Y a 2 +0.1 and that at the relevant column depth 
tr ^ tr 

t s y the dispersion is a 2 + 1/Y & a l + 0.05 (for Y £ 20). Since 

tr tr 'w 

a 2 « 1 this implies that most particles will have pitch angles less than 
tr 

It) degrees so that according to equation (3) most of the gamma-rays will 
also be direoted downward. Oaly a fraction of order of Y -2 of the 
gamr*. sprays will be in the visible hemisphere, and the Intensity of the 
feswsfc at the limb will be at most twice that of those occurring at the disk 
center. This cannot produce the observed distribution of figure (1) and 
will reduce the bremsstrahlung yield of gamma-rays by a factor of 

y~ z < 0.002. We can therefore rule out this highly beamed Injection also. 

The only remaining possibility is that a 2 r s a 2 B 0 /B tr so that at the 
transition region the electron distribution has become nearly isotropic in 
the downward direction, Belcw this region, because of this broad pitch 
angle distribution, analytic solutions are not possible (even though the 

field can be considered uniform and the synchrotron losses can be 

neglected) and a numerical treatment is necessary. However, we can 

estimate the amount of further dispersion in the diffusion approximation as 
follows. A distanoe Coulomb encounter produces an energy loss &y« 1 and 
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oau'»03 a deflection 5a = (25y/y 2 )^ . Thus a total of y/5y such 
encounters will produce a dispersion of the order of (5a) 2 = 2/y. For 
example, if at the transition region the electron distribution is isotropic 
(in the downward hemisphere), then without any dispersion the flares would 
be visible only within an angle 1/y of the limb, which as shown by the thin 
solid line in figure 1 is in disagreement with the observed distribution. 
However, with the addition of the dispersion caused by the Coulomb 

collisions, there will be some particles moving upward with pitch angles 
between ir/2 and ir/2 + (2/y) 5 * so that the expected distribution will be 

wider (heavy solid curve in figure 1) which is similar* to the observed one. 

We cannot be any more quantitative than this using analytic methods 

because, as described in the next section, the other factors mentioned in 
§1 must be taken into aocount before reaching any quantitative conclusions. 
However, to show the validity of the above arguments, in figure 2 we shew 
the expected directivity, at various relativistic photon energies, for 
Isotropic injection in a uniform field obtained from the exact 
(Fokker- Planck) numerical treatment of the electron transport and the exact 
evaluation of the bremsstrahlung emission (i. e. without the approximation 
in equation 2). This clearly demonstrates the decrease in the dispersion 
with higher energies. The portion of the 10 Msv curve for 0 > tt/ 2 is 

shown by the dashed line in figure 1 , which clearly is wider than the thin 

-h 

line and has a characteristic dispersion of ~ y . This is not in exact 
agreement with observation, indicating that the effects described below may 
be important. A more detailed analysis of the problem with the inclusion 
of the non-uniform fields and synchrotron losses will be described 


elsewhere 
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We conclude, therefore, that the observed distribution of the 
gamma-ray flares indicates that the pitch angle distribution of the 
electrons will be nearly isotropic in the downward direction and that the 
magnetic field convergence from the acceleration region (presumably in the 
corona) to the transition region is small. 


XV. THE INFLUENCE OF OTHER FACTORS 


Our discussion so far has been limited to the comparison of the 
directivity of the bremsstrahlung radiation, which we define as 


D(k,0) = J(k,0)/J(k) , J(k) 


U(k, 0)dJ) 


( 10 ) 


with the observed distribution (i. e. the number of gamma-ray flares 
brighter than some threshold flux) across the solar disk. Such a 
comparison is valid only approximately for some specific and limited 
intrinsic distribution of fluxes J(k) (i. e. for certain luminosity 
functions). For a more complete analysis of this problem a knowledge of 
the luminosity function is also required. 


In order to demonstrate this, let us assume that the luminosity 
function is described by a function i|/ (J) such that 


<P(J) = 


CO 

|" i|/(x)dx 
J 


( 11 ) 


is the rate of occurrence of gamma-ray flares with flux J(k) > j. Then the 
observed distribution of flares as a function of the heliocentric longitude 
H = ir - 0 in a sample limited to flares with fluxes greater than will be 

n. k (H) = <f>[J o /D(k,0)] 


( 12 ) 
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It is clear then that the distribution n will be proportional to 
D(k,0) qnly if 4>(J) * J -1 or ^(J) « J -2 . A large number of flares is 
required for an independent determination of the luminosity function. The 
present sample of gamma- ray flares is not sufficiently large for this task. 
However, at lower energies (photon energies of 25 keV), the luminosity 
function has the form $(X) « X -1 (cf. Lin 1984) so that if the same were 
true at 10 Mev the above comparison would be accurate. 


Another effect which must be taken into account is that the observed 
sample of > 10 Mev gamma-ray bursts has been selected from a larger sample 
of bursts which are selected for having fluxes greater than certain 
thresholds at lower (300 keV) photon energies. The simple comparison we 
have made would be valid if the distribution across the solar disk of this 
parent sample is uniform. Any non- uniformity in this distribution will 
affect the analysis carried out here. In fact, this distribution is shown 
to be non-uniform (Vestrand et al 1984) so that in absence of other effects 
the observed distribution at 10 Mev should be compared with 


n 


l o Mev 


(n) 


3 0 0 


Kev (n)<f)[ J o /D( 10 Mev, 0)3 


(13) 


As mentioned earlier, the obscured n (n) is much broader than the 

300 keV 

observed n _ (rj) 30 that the predominant factor in (13) is the 
l o Mev 

cumulative luminosity function $ at 10 Mev. 

A more complete analysis, which is beyond the scope of the present 
paper, will require knowledge of the directivities of bremsstrahlung at 
both 300 keV and at 10 Mev as well as the bivariate luminosity function 
ijj [J(300 keV), J(10 Mev)]. 


Page 


Finally, we should mention that any deviation of the magnetio field 
direction at large column depths (t ~ 20, N ** 10 24 cnf 1 2 ) from the local 
vertical will alter the transformatin between the angles r), 9, etc., and 
change the shape of the expected distribution. We should note, however, 
that the narrowness of the observed distribution indicates that the 
direction of the field lines cannot deviate strongly from the local 
vertical direction. 


V. SUMMARY AND CONCLUSIONS 

We have considered the observed distribution of flares with strong 
emission at photon energies of greater than 10 Mev in the framework of the 
model where the electrons responsible for the radiation are accelerated 
above the transition region and are Injected in a closed magnetic loop. We 
have shown that: 

1) Bremsstrahlung by relativistic electrons is the mechanism for 
production of these gamma-rays and that the production site must be deep in 
the photosphere where the predominant energy loss process is Coulomb 
collisions. The synchrotron losses can play a significant role if the 
magnetic field strength in the corona exceeds 1000 gauss, and the inverse 
Compton process could be important if the soft photon energy density in the 
flare region exceeds that of the quiet sun by three or four orders of 
magni tude. 

2) The magnetic field variation from the acceleration region to the 
photosphere must be small to prevent trapping of the particles in the 
corona where they will produce mainly microwave and millimeter wavelength 
photon via synchrotron radiation (and very little gamma rays) for a period 
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much longer than the observed duration of the gamma-ray bursts. 

3) The pitch angle distribution of the accelerated electrons cannot be 
beamed strongly along the field lines because this will result in the 
gamma-ray emission being directed toward the photosphere and, therefore, 
invisible from the Earth. 

The above three conclusions can be quantitatively expressed a3 
follows: If B q and B fcr are the strengths of the field lines at the 

acceleration region and below the transition region, respectively, then the 
characteristic value of th9 sine of the p.itoh angle of the accelerated 
particles must be of the order of (B /B ) • 11515 requirement of nearly 
uniform field and/or isotropic pitch angle distribution is similar to those 
deduced from the consideration of the spatial structure of the microwave 
radiation (Petrosian 1982). 

Finally, as shown in Section IV, a more quantitative conclusion in 
regard to the details of the model must await the determination of the 
distribution of the intrinsic fluxes of both gamma-ray and hard x-ray 
producing electrons. 


APPENDIX 

1) Synchrotron Losses: The radiation reaction force on a particle of 
mass m, charge e, energy yme 2 and velocity v = gc is (Panofsky and Phillips 
1962). 


F r = (2e 2 /3c 2 )[(yVg)g + y 2 g] 

~ *v »v ** rw 


(A.l) 


where n denotes time derivative. If the energy loss per gyro-period in a 
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magnetic field of strength B is small compared to the particle energy, then 
the particle will undergo helical motion with a pitch angle cosine y and 
with 

g * 3 = 0, 3 * -B(l-y 2 )(27rv b /Y) 2 * * , v fe = 2.8xl0 6 Hz(b/ 100 gauss). (A. 2) 


Til a rate of change of parallel and perpendicular momenta due to the 
radiation reaction then becomes 


P X 1 “ -(87r 2 e z v^/3c 2 )y 2 Bi 3 


11 * P, = -(8ire 2 v 2 /3c 2 )(Y 2 3l + 1)3. 


(A. 3) 


Transformation of these to rate of energy loss and pitch angle cosine 
change gives 

(Y/y) s = -(8Tr 2 r o v 2 /3c)Y3 2 (l-y 2 ), (y/y) s = (8Tr 2 r o v 2 /3c) (1 -y 2 )/ Y (A.4') 

These imply a continuous loss of energy and increase in the absolute 
value of y , At extreme relativistic energies the change in y is 
negligible while at non- relativistic energy the rate of change of energy is 
smaller than the rate of change of y , In general, $ = constant. 

2) Coulomb collision losses. The Coulomb losses are well known and 

for non- relativistic regime were given in LPI. A more complete description 

can be found in Leach (1984), where the various logarithmic factors are 

explicitly evaluated. Note that the extension of LPI result Into the 
relativistic regime gives an incorrect expression. For correct expression 
the quantity (3 +Y )/t in the coefficient C^ should be replaced by unity. 

In the present context we are interested in losses in a neutral gas. 
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Ignoring 20 to 40 percent contribution from helium, we can then write 

(y/y). =■ ~(4-rrr 2 c nlnA)/ey , (y/y)“ -(4Trr 2 c nlnA' )/0 3 y 2 (A. 5) 

where lnA= 21nA ; - 15 (for y « 20) and n is the neutral hydrogen 
density. For ionized hydrogen n would stand for the electron or proton 
density and InA - lnA x /2 « 20. Note that unlike energy, which in the cold 
ambient plasma always decreases, the pitch angle can either increase or 
decrease so that the higher order diffusion terms must be considered. The 
Fokker-Planok equation described below takes this into account. As evident 
at relativistic energies (y/y) c « (Y/Y) c • 

3) Inverse Compton Losses: The energy loss rate and the rate of 

change of pitch angle of relativistic electrons in a soft photon gas of 
energy density C is similar to that of synchrotron losses with the 
magnetic energy density B (1- y ) /8 tt replaced by C so that 

<Y/Y) x “ -(16tc 2 c/3) (Y£ Y /mc 2 ) , (y/y)^ - (y/y^/y 2 (A. 6) 

c 

As described in the text for the present application, the inverse Compton 
losses are negligible with respect to the synchrotron and Coulomb losses so 
that a more accurate analysis is not needed 

4) The Kinetic Equation: In a strong magnetic field the gyro- radius 

is much smaller than the mean free path and the scale of the spatial 
gradient of the field, plasma density, etc. Then the independent 
parameters of the distribution function are the time, distance of the 
particle guiding center along the field, s, and the two components 
(parallel and perpendicular) of the momenta. It is convenient to replace 
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the last two by the energy Es y - 1 and the cosine of the pitch angle, y. 
The Fokker- Planck equation for the distribution function of f(E,y,s,t) is 


3f 

c3t 


+ Sy|f-e 


dlnB 

2ds 


(i -y 2 ) 


3f 

ay 


3_ 

3E 


(Ef) 


_ a- 

ay 


(yf) 



(A. 7) 


In general, the higher order terms and the diffusion in energy can be 
neglected. The time dependence can be ignored if the acceleration and 
injection time scale of the electrons, as deduced from the modulation and 
duration of the observed emissions, is much longer than the microscopic 
time scale (such as energy loss time scale, etc.). Alternatively, for 
shorter events the integration of {A. 7) over the duration of the flare is 
equivalent to setting 3f/3t = 0 and having f = f(E,y,s) represent the time 
integrated distribution. The result of any emission calculation should 
then be compared with the time- integrated fluxes. 


The various possible solutions of this equation at non-relativistlc 
energies in uniform field and when the dominant 

process were discussed in LPI (see also Leach 1984). As described in the 
text, the observations dictate a flare plasma condition whereby the Coulomb 
collisions will be the dominant interaction of the accelerated electron. 
Furthermore, these interactions will take palce over a small distance at 
high photospheric densities where the field can be considered uniform. In 
that case equation (A.7), when written in terms of the flux F = c0f, 
reduces to (for steady state 3f/3t = 0) 





dT ~ 4irr 2 InAnds 
o 


(A. 8) 


I'tw - 


««. .u. ,Vl.» & 
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Integrating over all depths, we find 




pF Q (y,E) 


G(p,E) 


»cc 

F(E,p,x)dx 

0 


(A. 9) 


where F Q (E,y) is the distribution of the accelerated electrons injeated at 
T s 0 or s = 0. In the extreme relativistio case we can set 3 = 1 and E = y. 
Even with this simplification, convenient analytic solutions are 

possible only for small pitch angles: (l~y 2 ) = a 2 /2, 3/9p = (1/a)9/9a . 

Following LPI or Leach (1984), we find that (A. 8) gives the following 
simple solution for an assumed gaussian, strongly beamed distribution F q at 
injection. 


F Q (y,y) = 2 F q ( y) exp { -a 2 /a 2 }/a 2 , 

F(y,p,x) = 2 F q ( y+T ) exp {-a 2 /(a 2 +;)}/(oc 2 +?) , x, = y~ l x/(y + x) . 


(A. 10) 


As evident, the distribution remains gaussian at all depths but with 
increasingly wider dispersion. However, the dispersion is always small 
because £ < 1/y. If a 2 > y -1 , then for determination of 

G(y,vO it is easier to integrate (A. 10) approximately than solve the 
differential equation (A.9). For this we write a 2 + X = a 2 (1+£/ a 2 ) and 
expand (A.9), keeping only the first order terms in £/a 2 . Integration over 
all x then gives 


G(y,y) = 2G o (y+T)exp{-a 2 /(a 2 + £)} /(a 2 + £) 


(A. 11) 
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where 


G 0 (Y) 


oo 

F 0 (Y')dy' , X - Y _1 - <y~ 1 > 
Y 


and 


< y“ l > 


00 

F o (Y / )d Y / /Y ]/G o (y) . 
Y 


(A. 12) 


(A. 13) 


clearly £a1/y with the proportionately constant depending on the 
distribution F (y) . For example, if F o (y)=y“ n , then £ = 1/(ny) and 
G 0 (Y) 11 y" n+1 . 

For a 2 < y _1 the above procedure leads to a complicated integral 

0 'v 

which for a power law F q (y) « y n gives approximately 


G(Y ,M) « Y~ n+1 Uye" 0 ^] 


(A. 14) 


This is similar to ( A- • 1 1 ) when y is replaced by ct 2 + y o x with ? = 1/y. 
The important result from this is that for relativistic electrons the pitch 
angle dispersion of the integrated (or mean) distribution remains small if 
the dispersion at injection is small. 
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FIGURE CAPTIONS 

Figure 1 . The histogram shows the observed distribution of the flares with 
gamma-ray emission at >10 Mev across the solar disk; 8 s tt/2 - 0 where n 
is the heliocentric longitude (ri = 0 at the disk center, r| = tt/2 at the 
limb). 

The lower thin curve shows the approximate behavior of the 

bremsstrahlung cross section at relativistic energies with the small angle 
—2 

from 0+ 0 2 k 2 ) joined smoothly to the predominant term of the large angle 

2 

form, lT^sln 2 S/l6(1-cos0) ; k = 20 (10 Mev) and 6 is the angle between the 
electron and photon momenta. This curve can also be representing the 
variation of the flare brightness aaross the solar disk for bremsstrahlung 
from electrons with horizontal component of momenta. 

The heavy solid line is a gaussian fit with a dispersion in angle of 
“k -ls with k r 20. 

The dashed line is the calculated (numerical Fokker- Planck) variation 
of the gamma-ray flux across the solar disk at 10 Mev (from figure 2). The 
angular dispersion is ”k but the form is approximately exponential. 
The additional dispersion could come from other factors discussed in § IV. 

Figure 2, Directivity of hard x-ray and gamma-rays from exact numerical 
relations from a uniform loop with isotropic (downward direction) electron 
injection (cf, LPII for more detail) at three different values of photon 
energy k. 0 is the angle between the photon momentum vector and the 
downward vertical. Heliocentric longitude r) = it -0. 
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